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Abstract. We explore from a theoretical perspective the dynamical response of small water clusters, 
(H20) n H30 + with n = 1,2,3, to a short laser pulse for various frequencies, from infrared (IR) to ultra- 
violet (UV) and intensities (from 6 x 10 13 W/cm 2 to 5 x 10 14 W/cm 2 ). To that end, we use time-dependent 
local-density approximation for the electrons, coupled to molecular dynamics for the atomic cores (TDLDA- 
MD). The local-density approximation is augmented by a self-interaction correction (SIC) to allow for a 
correct description of electron emission. For IR frequencies, we see a direct coupling of the laser field to 
the very light H + ions in the clusters. Resonant coupling (in the UV) and/or higher intensities lead to 
fast ionization with subsequent Coulomb explosion. The stability against Coulomb pressure increases with 
system size. Excitation to lower ionization stages induced strong ionic vibrations. These maintain rather 
harmonic pattern in spite of the sizeable amplitudes (often 10% of the bond length). 
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1 Introduction 



The structure of water is generally considered to be rather 
ordered due to hydrogen bonds and with certain struc- 
tures prevailing. There is, in turn, a general agreement 
among models about the existence of small and medium- 
sized clusters of water molecules [T] . Furthermore, experi- 
mental evidence point to the fact that under appropriate 
conditions, water clusters may ionize. There thus exists, 
since rather long, a large body of investigations on charged 
water clusters, either anions [21314] . cations [5J or proto- 
nated clusters 17181. 



The dynamics of water clusters appears as a key issue 
in various scenarios, e.g., for the chemistry of formation of 
droplets and clouds in atmosphere [9]. Water clusters are 
also possible candidates as transient intermediate stages in 
liquid water [10] . They might thus play a key role in many 
chemical and physical processes |llll2ll3ll4U5j . The case 
of charged clusters, and especially ionized clusters, is par- 
ticularly relevant for studies on radiation damage, as wa- 
ter is the natural "environment" of biomolecules. This has 
triggered recent studies on the influence of a finite number 
of water molecules coating a biomolecule [15] . The direct 
analysis of the irradiation of protonated water clusters is 
also a key issue [T7] . 



From the theoretical point of view, water clusters have 
been considered since long at a structural level [21415) , for 
a recent summary see [18]. In recent years, studies of irra- 
diation of biomolecules, possibly in a water environment 
have also been started [19I20J , but usually without an ex- 
plicit time-dependent treatment of the ionization process 
itself. In the present paper, we aim at studying the re- 
sponse of small protonated water clusters to an irradia- 
tion by a short laser pulse. The choice of laser irradia- 
tion rather than excitation by a charged projectile is to 
some extent a matter of convenience as it allows for ex- 
ploratory studies a discrimination and a systematic scan 
of various dynamical regimes without dealing with compli- 
cations of scattering geometry. Furthermore, there exists a 
well documented literature on the laser irradiation of clus- 
ters, metallic or rare gas mostly [21122] in a wide range of 
excitations, which provides a robust background to study 
irradiation of other sorts of clusters such as water clus- 
ters. The paper is organized as follows. We first briefly 
recall the basic inputs of the model, in section [51 As a 
next step we first consider the optical response which is 
known to provide a key entry point to understand irradia- 
tion dynamics (section [3]). We then address irradiation by 
an intense laser filed, exploring in particular the influence 
of laser frequency and intensity as key characteristics of 
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the laser pulse (section^. We finally draw our conclusions 
in section [5] 



2 Model 

In order to perform microscopic simulations of dynam- 
ical processes, we employ time-dependent density func- 
tional theory (TDDFT) [23124125] for the valence elec- 
trons (described in terms of single-electron wave functions 
{ip a (r),a = l...JVei}) combined with classical molecu- 
lar dynamics (MD) for the ionic cores (described by their 
classical coordinates {Hi, I = 1 . . . -/V; on }). The ions are 
here H + , and 6+ . The starting point is an energy which 
includes the electron kinetic energy, the direct Coulomb 
energy, the exchange-correlation energy, the coupling en- 
ergy between electrons and ions (described in terms of 
non-local pseudo-potentials of the type and tuned 
to provide a most efficient description on the numerical 
grid |27j). and the energy coming from the interaction of 
the system with the external time-dependent laser field 
coupling to both, electrons and ions. The functional for 
the electronic exchange-correlation energy is taken at the 
level of the local-density approximation (LDA) with the 
actual form from [28] • The electronic energy is augmented 
with an average-density self-interaction correction (AD- 
SIC) US]- The ADSIC is a crucial ingredient to provide 
correct ionization potentials (IP) which, in turn, is neces- 
sary to describe electron emission correctly. 

Equations of motion are then derived in a standard 
manner by variation of the given energy. Variation with 
respect to tp a leads to the (time-dependent) Kohn-Sham 
equations for the single-particle wave functions. Variation 
with respect to the ionic coordinates yields Hamiltonian 
equations of motion for ions, thus treated by classical 
MD. Electronic and ionic equations of motion are solved 
simultaneously leading to coupled TDLDA-MD descrip- 
tion |3QI31j . The treatment goes thus beyond the adia- 
batic approximation. It becomes equivalent to the Born- 
Oppenheimer (adiabatic) approximation in regimes where 
the latter is valid. 

Time-dependent fields and wave functions are repre- 
sented on a 3D coordinate-space grid of dimensions 72 x 
72 x 64 or 72 3 , with a grid spacing of 0.412 a,Q. The elec- 
tronic wave functions are propagated in time by the time- 
splitting method [32]. The Coulomb problem is solved 
by a fast Fourier technique [33]. The ionic equations of 
motion are solved using the Verlet algorithm. Absorbing 
boundary conditions are used to properly remove outgo- 
ing electrons [3013 lj . Thus the total number of electrons 
N = N(t) decreases in time. The number of escaped elec- 
trons N csc = N(t = 0) — N(t) is a measure of average 
ionization. 

The electronic ground state wave functions are ob- 
tained by damped gradient iteration |34j . Ionic ground 
state configurations are obtained by cooling pseudo-dyna- 
mics [3S] . Results for the structure of a few selected water 
clusters are shown in figure [T] Gray balls stand for the 
oxygen atoms and are labelled 01, 02, etc. White ball, 
labelled HI, etc., represent the hydrogen atoms. The O-H 
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Fig. 1. Schematic view of the initial configuration of the sys- 
tems H2OH3O+ (H 2 0) 2 H 3 + , and (H 2 0) 3 H 3 + . Note that 
all indicated frames are positively oriented. 



bond lengths are also indicated and take basically three 
values, around 1.8, 2.0 and 2.7 oq. 



3 Optical response 

The spectral distribution of optical absorption strength 
(optical response) provides key information on the cou- 
pling of a system to (laser) light. The initial response is 
mostly mediated by electrons. The optical response then is 
related to eigenfrequencies of electronic excitations which 
exhibit well isolated peaks at frequencies below emission 
threshold and continuum above. These eigenmodes are im- 
portant because they represent the doorways for laser ex- 
citation. In the case of metal clusters, the major absorp- 
tion strength lies in the optical domain with a marked 
collective character. They are known as the Mie plasmon 
resonances and correspond to a collective oscillation of the 
electron cloud with respect to the ionic background, with 
rather well established systematic behaviors with cluster 
size [31l36j . The situation is more involved in organic sys- 
tems where varying bonding can lead to large differences 
from one cluster to the next. Furthermore the typical 
eigenfrequencies lie in the UV part of spectrum rather 
than in the optical domain and thus require laser pulses 
with higher frequencies (single photon) or higher intensi- 
ties (multiphoton processes). 

The optical absorption strength can be computed eas- 
ily within our TDLDA description by the technique of 
spectral analysis 37 38 , for a detailed discussion see |39j . 
To that end, we apply initially an infinitely short pulse 
covering all frequencies with equal weight, propagate the 
system for a sufficiently long time (depending on the de- 
sired spectral resolution) , record the time-dependent dipole 
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momentum D(i), and finally Fourier transform the sig- 
nal into the frequency domain. The strength is then ob- 
tained as the imaginary part of D(w). Results for opti- 
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Fig. 2. Spectral distribution of optical absorption strength for 
the systems H 2 0, H 3 + , H 2 OH 3 + , and (H 2 0) 2 H 3 + . The 
corresponding IP are indicated by an arrow. 



cal absorption spectra of H2O, H3O" 1 ", H20H 3 + , and 
(H20)2 HaO + are shown in Figure [2j The spectra exhibit 
typical structures with well marked and isolated eigenfrc- 
quencies at "low" frequencies below the ionization poten- 
tial (IP) for the given system and a trend to a more fuzzy 
" continuum" response at higher frequencies above the IP. 
The figure shows that infrared and visible light is below 
any significant peak and thus off resonant. Light pulses in 
the UV range can couple more or less strongly, dependen- 
ing an the detailed frequency. The figure does also indi- 
cate the IP for each system. These are in detail: IP(H20) 
= -14.4 eV, IP(H 3 0+) = -22.3 eV, IP(H 2 OH 3 0+) = 
-20.6 eV, and IP((H 2 0)2H 3 + ) = -16.9 eV. 



4 Irradiation by a strong laser field 

We now turn to direct analysis of the irradiation and re- 
sponse of small protonated water clusters to short and in- 
tense laser pulses. Even in the UV regime, the laser wave- 
length remains much larger than the size of the system 
under consideration. So the long wavelength limit applies 
such that the dipole operator describes the spatial part 
of laser pulse. Thus we model the laser pulse as the time 
dependent external field 

Mas = E e poi ■ r cos(ut)f(t) (1) 

with the amplitude E$ being related to the intensity I as 

E (Ry.ao x ) = 1.07 x 10- 8 -y//(W/cm 2 ) and (linear) polar- 
ization e po i. In all our calculations, we use e po i = e x (see 
Figured]). We will see below, in the discussion of figure 
|4l that the laser polarization has only a limited impor- 
tance for the basic mechanisms. The laser frequency u> is 
carried in the cos term. The pulse profile f(t) is chosen 



as a sin 2 characterized by its Full Width at Half Maxi- 
mum (FWHM). In all calculations below we have chosen 
a FWHM of 20 fs and we have focused on the simple case 
of (H20)H30 + , exploring only at some places larger clus- 
ters. 



4.1 Dependence on frequency 
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Fig. 3. Irradiation of (H 2 0)H30 + by lasers with polarization 
along x axis, FWHM of 20 fs, intensities of about lo = 10 13 
W/cm 2 , and three frequencies as indicated. The number of es- 
caped electrons N csc is also reported for each case. As a func- 
tion of time : Electronic dipole moment in x direction (left 
column), and ionic x coordinates (right column). The thick 
curves stand for the position of the oxygen atoms, while the 
thin lines show the x coordinates of the hygrogen atoms. 



Wc first consider the influence of laser frequency on 
irradiation dynamics of (H20)H30 + , exploring a varia- 
tion of laser frequency at basically constant laser inten- 
sity. Figure [3] shows results for three representative cases: 
a frequency in the infra-red (IR), another one in the visi- 
ble range, and an ultra-violet (UV) one. The IR and vis- 
ible frequencies lie below the range of typical electronic 
eigenfrequencies in (H 2 0)H 3 + (see Figure [5]), while the 
UV case lies in the middle of the eigenfrequencies. The 
overall dynamical behavior is strongly related to the laser 
frequency. We observe a typical out-of-resonance behav- 
ior for IR and visible frequencies where the amplitude of 
the dipole signal (left column of Figure |3]) is proportional 
to the envelope of the laser field. There is practically no 
electron emission for the given laser intensity. The UV 
frequency, being in resonance with some of the systems 
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excitation modes, displays a much different behavior: the 
strong dipole response lasts longer than the laser signal 
with a sizeable ionization. Correspondingly, the overall 
behavior of the cluster dynamics behaves quite differently. 
The right column of Figure [3] displays the time evolution of 
the ionic positions of constituting ions along the laser po- 
larization axis. While the two off-resonance cases exhibit 
only low amplitude ionic vibrations, the UV case leads to 
an almost immediate explosion of the (H20)H30 + . 

Another interesting aspect concerns the difference be- 
tween the visible and IR cases. In the frequency in the vis- 
ible range, one can hardly see any ionic perturbation while 
the IR irradiation leads to clear ionic vibrations with an 
small, but not vanishing, amplitude. This difference can 
be traced back to a direct coupling of the ions (predomi- 
nantly the hydrogen ions) with the laser itself. This seems 
a bit surprising in view of the rather small FWHM of the 
laser pulse (20 fs). But the small hydrogen mass allow 
such a rapid coupling. This coupling is observed in the IR 
domain and not for the visible light because typical ionic 
vibration frequencies in small water clusters precisely lie 
in the considered IR frequency domain [40) . Estimating 
the H vibration frequency in this case leads to a value of 
hu ~ 0.41 eV, in fair agreement with previously published 
results which shows peaks around 0.47 eV [ID]. Closer in- 
spection shows that the IR laser pulse simply excites inter- 
nal vibrations of hydrogen atoms mostly along OH bonds. 

In order to explore this effect in more detail, we now 
focus on the distances between O and H atoms. Indeed, 
because of their larger masses, O atoms basically remain 
fixed during the short time scales explored here and we 
can thus focus the further analysis on O-H bonds only. As 
can be seen from Figure [1] there arc basically two types of 
O-H bonds, the ones between O and external H and the 
internal ones along the sequence H-O-H. The associated 
bond lengths are different and so are the corresponding vi- 
bration frequencies. It is thus interesting to consider these 
two classes of cases separately. The results are plotted in 
Figure 01 One should first note the remarkable regular- 
ity of oscillations, especially for external H atoms (upper 
panel), with about all having the same frequency. One 
also clearly observes oscillations in opposite phases de- 
pending on the geometry. This is particularly obvious for 
the two internal bonds (lower panel of Figure U|) . Different 
heights of amplitudes indicates a different strength of the 
response which is predominantly a geometry effect due to 
the orientation of the bond relative to laser polarization, 
here the x direction. Indeed, the Oi-Hi and O2-H5 bonds 
are the external O-H bonds which are aligned with the x 
axis the most (see Figure [I} and which exhibit the largest 
amplitudes of oscillations. Changing the laser polarization 
(along the y or the z direction for instance) would produce 
a largest excitation of other O-H bonds but the mecha- 
nisms of excitation would remain the same. Finally one 
should note the amplitude of oscillations which may be 
rather large, typically 10% of the actual bond length but 
without destroying the harmonicity of the oscillations, the 
latter aspect being especially true for external H atoms. 
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Fig. 4. Irradiation of (H20)H30 + by a laser with polarization 
along x axis, FWHM of 20 fs, and intensity and frequency as 
indicated. The number of escaped electrons is also given. Are 
plotted as a function of time the O-H bond lengths : the short- 
est ones (top), according to the atom numbering of Figure [TJ 
and the longest ones (bottom). 



4.2 Dependence on intensity 



We have identified in Figure|3]the resonant coupling of the 
laser pulse with internal IR vibrations along OH bonds 
and explored it in more detail in Figure |U Figure [5] ex- 
plores how this behavior depends on laser intensity. It 
shows the time evolution of the distances between the 2 
oxygen atoms which lies around 4 ao and the typical OH 
bonds whose length is about 2 ao. The influence of laser 
intensity on the ionic response is obvious. With increas- 
ing laser intensity, we observe increasing ionization. Up to 
about net ionization one, there is little effect on the general 
vibrational pattern. They gain, of course, amplitude due 
to the larger field strength. But on the time scales com- 
puted here, the systems seems to resist fragmentation and 
keep its vibrational response. Only for ionization substan- 
tially larger than one, the dynamical evolution changes 
qualitatively to the extent that the now large repulsion 
leads to an immediate fragmentation of the cluster and, 
of course, disappearance of the vibrational response. All 
in all, we thus observe a remarkable robustness of the vi- 
brational response with increasing laser intensity. 
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Fig. 5. Irradiation of (H20)HaO + by lasers with polarization 
along x axis, FWHM of 20 fs, frequency of 0.48 eV, and six dif- 
ferent intensities varying about Io = 10 14 W/cm 2 as indicated. 
Intensities increase in clockwise order. The average ionization 
(number of escaped electrons N eac ) is also indicated for each 
case. Various bond lengths (0-0 and shortest O-H) are de- 
picted as a function of time. 



4.3 Dependence on size 

It is finally interesting to explore how dynamical evolu- 
tion depends on cluster size. We consider here the case of 
(H20) n H 3 + clusters for n = 1,2,3. We again take the 
same IR frequency at 0.48 eV and a moderate laser inten- 
sity to remain in the vibrational regime in (H20)H30 + . 
Various OH bond lengths are plotted as a function of time 
in Figure O Because of the quickly increasing number of 
OH bonds with increasing size, which makes a figure in- 
cluding all of them unreadable, we have selected a few rep- 
resentative cases, following the labelling given in Figure [TJ 
Note that, although we use the same laser in all cases, the 
net ionization differs from one case to the next, due to 
different ionization potentials: —20.6 eV, — 16.9eV, and 
— 15.7 eV for n=l, 2, and 3 respectively; associated net 
ionizations are 0.34, 1.44, 2.00 respectively for n = 1, 2, 3. 
This should not be a problem as we have seen that the 
effects are to a large extent robust with net ionization. 
Indeed, at least on the short time scale we have consid- 
ered here, we observe in all cases comparable vibrational 
response, in spite of the large ionization in the largest clus- 
ter. Note, however, that the larger ionization appears in 
the larger system which has naturally a higher capabil- 
ity to cope with a high charge state. Such an effect has 
been found also for the Coulomb fragmentation of metal 
clusters [4"T] . 

Figure O exhibits interesting features in detail. The 
upper panel compares typical external O-H vibrations in 
the n = 2 and n = 3. In both clusters, the corresponding 
bonds have very similar lengths and frequencies, as is clear 
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Fig. 6. Irradiation of (H20)„H30 + , for n — 2, 3 by lasers with 
polarization along x axis, FWHM of 20 fs, frequency of 0.48 eV, 
and intensity of 2.6 x 10 14 W/cm 2 . Various O-H bond lengths 
are depicted as a function of time, with labels according to the 
atom numbering of Figure [1] 



form the plot. The responses are remarkably similar. The 
middle panel focuses on (H20)2H30 + internal H-O-H se- 
quence comparing H-0 and O-H vibrations, which turn 
out to be pretty similar and nicely in phase opposition, 
simply because the O atoms basically do not move. Finally 
the lower panel concerns (H20)3H30 + , again focusing on 
a H-O-H sequence and with the same conclusions as in 
the middle panel. 



5 Conclusion 

We have explored in this paper the dynamical response 
of small water clusters (in particular (H 2 0)H 3 + ) to a 
short and intense laser pulse for various frequencies and 
intensities covering the regime of stable vibrations up to 
Coulomb break-up of a highly ionized cluster. The first 
stages of response of the system are, as is well known from 
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other systems, primarily taken up by the electrons and we 
have recovered this behavior. For low frequencies in the 
infra-red (IR), we also see a direct coupling of the laser 
pulse to the ionic centers of the clusters, particularly to the 
most mobile H + ions. This coupling becomes sizeable be- 
cause the ionic vibrations are close to resonant conditions 
for IR frequencies. The effect is robust with respect to 
change of laser intensity, even for sizeable extra ionization 
up to typically one more charge unit. We have also stud- 
ied larger hydrogen clusters with one or two more water 
molecules and found similar behaviors. We see here that 
the limits of stability against Coulomb break-up is shifted 
to higher ionization with increasing system size. The state- 
ments about stability have to be taken with care because 
our observation times are rather short. This means that 
we have discussed at least a stable transient stage. One 
needs to follow the evolution for longer time spans to find 
out more about the long-time stability. Work along that 
line is in progress. 

This work was supported by the Deutsche Forschungsgemein- 
schaft (RE 322/10, RO 293/27), Fonds der Chemischen Indus- 
trie (Germany), a Bessel- Humboldt prize, a Gay-Lussac prize, 
the French Agence Nationale de la Recherche (ANR-06-BLAN- 
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